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8 Field of the Invention 

9 This invention generally relates to a device and method for monitoring fluid properties, 

10 and specifically relates to a density/solute monitor having ultrasound probes for continuous 

1 1 monitoring of the ultrasound velocity of fluid in a biological or chemical processing system in 

12 order to determine fluid density, compressibility, solute concentration, and the fluid flow and a 

13 method for using the same wherein the probes are integrated with modality measurements such 

14 as optical absorbance, conductivity, impedance, magnetic resonance, radiation attenuation, and 

15 tracers of fluid. 
16 

1 7 Background 

1 8 Two classical methods of measuring density of a fluid include: 1) measuring the weight 

1 9 of fluid in a flask of fixed volume; and 2) employing the buoyancy of a density float for the 

20 assessment of fluid density. Both of these methods require collection of large samples from a 

21 fluid-processing device such as a pipeline or reactor for off-line measurements. For a given 

22 solution, density relates to the solute concentration of the solution. Although the density 

23 measurement is not specific to what solute is in the solution, these two density measurement 
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1 methodologies and others to be described later have been used as a means to assess solute 

2 concentration. 

3 A mechanical device based on resonance has been available to measure density of a fluid 

4 sample or that of a flowing fluid on-line. The device has a hollow U-tube with its two ends 

5 fixed on a heavy base. The fluid can be infused to fill the U-tube or made to flow along the tube. 

6 By measuring the frequency that the U-tube resonates, one then determines the mass of fluid in 

7 the U-tube. Since its volume is fixed, the mass is converted to the fluid density. This 

8 mechanical density measuring system (MDMS) has high sensitivity and reproducibility in the 

9 dynamic measurement of fluid density. 

10 As a fourth density measurement method, one measures the sound velocity of fluid for 

1 1 the determination of the compressibility and density of fluid. Krivitski, in U.S. Patent Numbers 

12 5,453,576 and 5,685,989 describes an apparatus and method for measuring several hemodynamic 

13 parameters by using a sound velocity sensor. The ultrasound transducer is excited to emit a pulse 

14 of ultrasound. After its passage through a fluid medium such as the blood, a receiving transducer 

1 5 senses the ultrasound pulse. A protocol to compare the excitation and receiving ultrasound 

16 signals determines the transmission time through the blood and subsequently its sound velocity. 

17 The information contained in the '576 and the *989 patent is incorporated by reference as though 

1 8 cited in its entirety. When their device is used to measure blood density change for the 

19 computation of blood volume, a linear approximation of a non-linear relationship is employed to 

20 convert the sound velocity to the density of the blood. Furthermore, the device of the *576 patent 

21 has limited sensitivity so as to require the imposition of a large change in blood density for 

22 accurate assessment of blood volume. 
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1 The system patented by Schneditz in US Patent 5,830,365 also utilizes sound velocity for 

2 the measurement of transmission time delay through the blood and then its total protein 

3 concentration. A clinical protocol to change the ultrafiltration rate as a patient undergoing 

4 hemodialysis treatment is described to produce the change in density, which is assessed through 

5 a sound velocity monitor. An equation is deduced to compute from the change measurement the 

6 blood volume circulating in the patient. The monitor to measure density is about one order of 

7 magnitude less sensitive than that provided by the MDMS or our density/solute monitor. As a 

8 result, the application of Schneditz' s method to measure blood volume is limited to cases that the 

9 change in sound velocity being imposed through the clinical protocol is large. 

10 The fifth density measurement method employs the attenuation due to the absorption of 

1 1 radioactivity by the fluid as a means to assess its density. Approval by regulatory agency is 

1 2 required for this method. 

1 3 Summary of the Invention 

14 The present invention is directed to a density/solute monitor including an ultrasound 

15 probe and a signal processing unit for accurately and reliably determining the phase shift of 

16 ultrasound transmission through a fluid and then the sound velocity of the fluid and a method of 

17 applying the same to biological or chemical processing systems. A set of equations and 

1 8 measurements by other modalities are incorporated to deduce from the sound velocity the 

19 compressibility, density, concentration of specific solute, and concentration of particulate matters 

20 of the fluid. The monitor can be used to determine the mass flow of solute, to improve the 

21 performance of chemical processes, and to optimize process design. The improvement and 

22 optimization can lead to more efficient collection of solutes, more solute purity in the collection, 

23 and better efficiency of the chemical processes. 
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1 The novel embodiments of the density/solute monitoring system include a signal 

2 processing unit with simple hardware and software to determine at high accuracy phase shift and 

3 transmission time of ultrasound signals; a two-fluid calibration procedure to convert the phase 

4 shift and transmission time in terms of sound velocity; an appropriate placement of the probe to 

5 time the passage of certain solute injected upstream of the probe; the use of two probes in series 

6 to measure volume flow; a procedure to work with the MDMS for the establishment of an 

7 empirical relation between the density and sound velocity of fluid and to account for the 

8 dependence on temperature; the calculation of the compressibility of the fluid to derive its 

9 relation with the sound velocity and density of the fluid; a set of computer files and equations 

1 0 specific to given solute, solution and density/solute monitor on the conversion of density to the 

1 1 solute concentration in the solution. 

1 2 The ultrasound probe can work alone or in combination with other detection modalities to 

13 achieve more functionality for the density/solute monitor. Other detection modalities include: 

14 (1) The use of optical absorbance and/or reflectance of light at frequency ranging from 

15 infrared to ultraviolet, impedance and conductivity of microwave, and absorbance of 

16 radiation for better identification of the solute of interest; 

17 (2) A procedure using the injection of certain solution and the sensitivity of the 

1 8 ultrasound probe to calibrate and determine the sensitivity of other detection 

19 modalities; 

20 (3) A procedure to detect the movement of tracers for system characterization. 

21 The use of these embodiments will enable one to achieve at least one of the following features: 

22 (1) A multi-functional detection system of low cost; 

23 (2) Accurate assessments on the flow and passage of specific solute; 
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1 (3) Efficiency in the collection and purification of specific solute with the technology of 

2 chromatography; 

3 (4) Crucial information for industrial engineers to optimize the process design; 

4 (5) Information for diagnosis and prevention of cardiac deficiencies in patients. 

5 By making use of the high sensitivity of the ultrasound probe, the multi-modality monitor 

6 gains additional power to determine the concentration of specific solute in the solution, the 

7 passage of specific solute through a chromatography column, and the dynamic changes of the 

8 solute in chemical or biological processing systems. In applying these probes and methods to 

9 paper, petroleum, chemical, pharmaceutical, food and bioprocessing industry; the engineers can 

10 determine more accurately the mass flow being transported through pipeline, achieve better 

1 1 solute purity in solution collected from chromatography column, and control more responsively 

12 the chemical or biological processing. The multi-modality methodology is applicable to human 

13 for determinations of blood parameters, vascular functions, and cardiac performance. The 

14 information provides key measurements for physicians to maintain homeostasis of the patient 

15 and to diagnose or to prevent cardiac deficiencies such as hypotension and shock in patients 

16 undergoing hemodialysis treatment or subject to trauma or burns. 

1 7 Multi-modality probes and methods are described: 

18 1 . To monitor the phase shift between the emitting and receiving ultrasound and the 

19 time of sound transmission in the fluid in pipelines or reactors; 

20 2. To use a set of relationships and procedures to convert these phase and time 

21 measurements into accurate assessment of sound velocity, density, compressibility, 

22 solute concentration, and their changes; 
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1 3. To monitor a number of fluid properties in optical absorbance, reflectance, 

2 conductivity, impedance, magnetic resonance, radioactivity attenuation, and tracers to 

3 better specify the solute being assessed and to better time the flow and passage of the 

4 solute; 

5 4. To employ these assessments for more efficient operation and control of chemical, 

6 physical and biological process common to chemical, pharmaceutical, food product, 

7 paper, and petroleum industries; 

8 5. To improve the probes and algorithm for use as blood volume monitor claimed in U.S 

9 [[p]]Patent 1 0/271 >686 numbers 6.740,036 and 6,485,427 B 1 . 

1 0 Brief Description of the Drawings 

1 1 Figure 1 is a schematic drawing of an embodiment of the ultrasound probe in the 

1 2 insertion mode. It has a pair of ultrasound transducers and a thermistor. The probe is inserted 

1 3 into a pipeline or reactor for measurements of fluid inside. In this design, the transducers and 

14 thermistor are in direct contact with the fluid. 

15 Figure 2 is a schematic drawing of an embodiment of the ultrasound probe in the clip-on 

1 6 mode. It has a pair of ultrasound transducers and a thermistor. The probe is clipped on to a tube 

17 or column with flowing solution. In this design, the transducers and thermistor have no direct 

1 8 contact with the solution. 

19 Figure 3 is a schematic drawing of an embodiment of the ultrasound probe in the cuvette 

20 mode. Pair of ultrasound transducers and a thermistor are mounted on the sides of the cuvette for 

21 which it contains the fluid for measurement. In this design, the transducers and thermistor are in 

22 direct contact with the solution. The housing containing the cuvette has the function of 

23 maintaining temperature at a preset value. 
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1 Figure 4 is an illustration identifying the phase and time shift among the excitation signal, 

2 that to the emitting ultrasound transducer, and that from the receiving ultrasound transducer. The 

3 excitation depicted has a constant frequency. Only over a finite time (five oscillation periods are 

4 depicted here and then repeated 26 cycles later), the excitation signal is passed through to excite 

5 the emitting transducer to emit ultrasound. After the transmission of ultrasound through the 

6 fluid, the receiving transducer picks up the ultrasound. Its amplified version is depicted here. The 

7 transmission time, exemplified by the time for valley E in the emitting signal to be transmitted to 

8 valley R in the receiving signal, is composed of n oscillation periods (an n of 15 is depicted) and 

9 a phase shift between the excitation signal and receiving signal ((p). 

10 Fig. 5 depicts the hardware employed to digitize the excitation signal (or the signal to the 

1 1 emitting transducer) and the signal from the receiving transducer, the storage of the data in the 

1 2 memory of interface processor, and then its transfer to the computer for the determination of the 

1 3 phase shift. 

14 Figure 6 is the density of saline, phase shift, and salt concentration in saline at 23°C, The 

15 density is measured by a MDMS and the phase shift by an insertion ultrasound probe. No gating 

16 of the excitation signal was employed for this data set. The range of density is achieved by 

1 7 varying the salt concentration in the saline. 
18 

1 9 Description of the Preferred Embodiments 

20 Figure 1 shows an embodiment of the ultrasound probe in the insertion mode. In this 

21 embodiment, the probe 21 includes an inserting mechanism 26, which is mounted with the 

22 ultrasound emitting transducer 22, the receiving transducer 24, and the thermistor 66. With the 
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1 insertion of the probe into a fluid system such as a chemical system or pipeline, the surfaces of 

2 the transducers 22 and 24 and thermistor 66 are directly in contact with the fluid flowing there. 

3 Transducers 22 and 24 are part of said ultrasound probe, which is attached to a signal 

4 processing unit. The signal processing unit is comprised of a function generator 20, an amplifier 

5 68, a dual channel analog-to-digital (A/D) converter 70, an interface processor 72, and a 

6 computing mechanism 74 as shown in Figure 5. The function generator 20 transmits a power 

7 signal, preferably via a cable, to activate the emitting transducer 22 into producing a train of 

8 ultrasound wave at an appropriate frequency (fuitrasound). One preferred frequency is about 5 

9 Megahertz, but any frequency deemed appropriate by one skilled in the art would suffice. 

10 This power signal is also digitized as an excitation signal by one channel of the dual 

1 1 channel (A/D) converter 70 within the ultrasound signal-processing unit. When the excitation 

12 frequency of the excitation signal is chosen as 5 Megahertz, the emitting transducer emits 

13 ultrasound at the frequency of 5 Megahertz. The sampling frequency of the A/D converter is 

14 chosen to be an integer multiple of the excitation frequency. This integer is designated as m. 

15 For an ultrasound/excitation frequency of 5 Megahertz and an A/D converter's sampling 

16 frequency (f sampling) of 65 mega samples per second (MSPS), m is 13, meaning there are 1 3 

17 digitized samples over one full ultrasound oscillation. Or, if an A/D converter with a sampling 

1 8 frequency of 105 MSPS is used, m will be 21 . When low cost A/D converters at higher sampling 

19 frequencies become available, the ultrasound frequency or the number of samples per oscillation 

20 cycle can be increased to improve the resolution of the density/solute monitor. The receiving 

21 transducer 24 receives the ultrasound wave after its passage through the fluid. The signal is 

22 amplified by amplifier 68 and digitized through the other channel of the dual channel A/D 
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1 converter 70. Both the digitized excitation and receiving signals are sent to the interface 

2 processor 72 for storage and a computing system 74 for analysis. 

3 The system of the present invention uses continuous measurements of phase shift and 

4 transmission time of ultrasound waves in a fluid to assess sound velocity, density, solute 

5 concentration, compressibility, and changes in these quantities. Once the ultrasound probe is 

6 inserted into a fluid and the signal processing unit described above transmits and digitizes 

7 ultrasound waves through the fluid, the computing system 74 is ready to determine phase shift 

8 and transmission time. The transmission time (Tshift) depicted in Fig. 4 is the time for valley E of 

9 the emitting signal after transmission in becoming valley R of the receiving signal. With the use 

10 of continuous wave, the period of each oscillation (To) is 1/ fuitrasound- To facilitate the 

1 1 determination of transmission time, it is divided into two parts. The first part is composed of n 

12 periods of oscillation, which is the time for valley E to reach a valley of the excitation signal that 

1 3 is closest to the valley R. A peak detector and a clock determine this shift. The second part 

14 reflects the phase shift from the valley of the excitation signal to the valley R. To be determined 

1 5 from the procedure described next, this phase shift [[f]]SL has the range in between -1 80° and 

1 6 1 80°. Accordingly the transmission time is expressed as: 

1 7 Tshift = To (n + [[f]]<a/360°) Equation 1 

1 8 or = (n + [[f]](a/360°)/fuitrasound Equation 2 

19 Let us identify the excitation and receiving samples stored on the interface processor as Ei 

20 and Ri respectively with i being the sampling index. Once some 1000 to 2000 samples from 

21 each signal are stored, the computer instructs the interface processor to transfer the data for 

22 processing. First, their average is determined and subtracted to achieve a zero average. Then the 

23 data are multiplied and summed as specified by Equation 3 and 4. 
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Ml = [[S]]S(Ei Ri) with the summation from i = 1 to N 



Equation 3 
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M2 = [[S]]2(Ei Ri+j) with the same summation as above 



Equation 4 



3 where j is chosen so that the receiving signal is shifted by about one quarter of a cycle. It is 



5 The function Round stands for the nearest round off of a number to an integer. The total number 

6 N used for the summation is chosen to be an integer multiple of m and to cover most of the 

7 period over which we have values for the receiving signal. When the number of cycles to be 

8 covered is larger than 60, our computation results indicate good sensitivity in sound velocity 

9 determination will be achieved. With Mi and M2 so calculated, we determine the phase shift of 

10 the receiving signal from the emitting signal ([[f]](a) as: 

1 1 [[f]](a = tan * [M2/(Misine) - cotO] Equation 6 

12 where [[f]](g and 0 are expressed in the unit of degree and 0 is 360° • (j • fuitrasound/fsampiing). If rn 

13 is a multiple integer of 4, then 9 = 90"" and Equation 6 reduces to the one commonly used in 

1 4 phase lock computation : 

15 [[fl](£ = tan'^[M2/Mi] Equation 7 

16 Two fluids typical to certain fluid processing will be used to calibrate the probe in the 

17 factory or in situ. Let the sound velocity of the two calibrating fluids be Ci and C2 and the 

1 8 corresponding phase shift be [[f]](gi and [[f]](E2. The time for the ultrasound to transmit from the 

19 emitting transducer through the fluid to the receiving transducer relates the measured phase shift 

20 by these equations: 

21 (n + [[fl](Ei/360°) To = L/ci Equation 8 

22 (n + [[f]]ce2/360°) To = L/C2 Equation 9 
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j - Round(m/4) 



Equation 5 
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1 where L is the distance between the transducers for this insertion probe. Suppose the sound 

2 velocity of the fluid designated for the measurennent is C3 and the measured phase shift is [[f|]c£3. 

3 Then they are related by 

4 (n + [[f]]ce3/360^) To = L/C3 Equation 1 0 

5 Equations 8, 9 and 10 can be reorganized to yield Equation 1 1 to determine C3 from [[f]]<a3: 

6 C,/C3 = 1 - (1 - ci /C2)([[f|]c£3 ' [[f]]cei)/([[fl]ce2 " [[f]]cei) Equation 1 1 

7 Using a series of fluid samples having a range of solute concentration, we can use the 

8 ultrasound probe and the MDMS to determine the sound velocity (c) and the density ([[?]]fi) 

9 respectively. It is known that the compressibility ([[?]]k) of the fluid relates to sound velocity 

1 0 and density by 

1 1 [[?]]k = [[?]]fi/c^ Equation 1 2 

12 By plotting the measurements and calculations against each other, we obtain a set of empirical 

1 3 equations for converting the measured sound velocity in terms of density, compressibility, or 

1 4 solute concentration. 

15 For most cases in industrial and clinical application, the difference among the three sound 

16 velocities is smaller than a few percentages. Thus, we can linearize Equation 1 1 to relate the 

17 sound velocity to phase shift by: 

18 C3= c, + (C2 - ci)([[f]](£3 - [[f]]cai)/([[f|]ca2 " [[fllseO Equation 13 

19 Since the change of density is also smaller than a few percentages, the change in density and that 

20 in sound velocity can be related through a linear equation. Accordingly, Equation 13 can be 

21 converted to the following form for the determination of density: 

22 [[?]]C3= [[?]]& + ([[?]]fi2 - [[?]]ei)([[f]]ffi3 - [[f]](ai)/([m]ffi2 - [[f]](Bi) Equation 14 
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1 where [[?]]fi3 is the density being measured, and [[?]]ei and [[?]]e2 the density of the calibrating 

2 fluids. In the case of protein solution, the density is linearly related to the concentration of 

3 protein C. If the fluid in the pipeline or chemical reactor also has its solute concentration linearly 

4 related to the density, we have Equation 1 5 to derive from the phase shift measurements the 

5 solute concentration C3: 

6 C3=Ci + (C2-Ci)([[fl]c£3- [[f]](Bi)/([[fl](a2 - [[f]]cai) Equation 15 

7 In the special case that the solute concentration for one calibrating fluid Ci is zero. Equation 15 is 

8 simplified to: 

9 C3 = C2 ([[f]]ffi3 - [[fl](&i)/([[f]]^ - [[f]]cei) Equation 16 

1 0 Equation 1 3 is regarded as a two-constant calibration equation of the monitor probe or 

1 1 probes of the monitoring system to convert phase [[f]](a3 to C3. In this equation, ci is given and 

12 the two constants are [[f]](ei and (C2 - Ci)/([[f]](e2 - [[f]]s&i). These two constants are determined 

13 by the two fluid calibration procedure. In the same way, one can define the two constants in 

14 equation 14 or 15 for converting phase to density or solute concentration respectively. 

1 5 Similar A signal processing procedure similar to the probe calibration discussed above is 

16 applicable to the case [[that ]] where the emitting ultrasound is identical to the excitation signal, 

17 i.e. without the gating shown in Fig. 4. In this scheme, wave reflection will take place and the 

1 8 total transmission time could not cannot be determined. However, we can still use Equations 3, 

19 4 and 6 can be used to determine the phase shift of the receiving ultrasound from the emitting 

20 one. With the values of [[f]](Bi and [[f]](B2 determined for a given ultrasound probe in the 

21 factory, one finally uses Equation 14 to 16 to determine from the phase shift measurement [[f]](e3 

22 the sound velocity, density, and solute concentration. This computation scheme as applied to 

23 blood, saline and plasma has been described in U.S. Patent 10/27^,0 8 6 filed October 18, 2002 
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1 No. 6,74Q.036 > to which the present application claims priority. The computation scheme and 

2 monitoring system can be used to determine sound velocity, density, and solute concentration 

3 measurements for various fluids in industries including, but not limited to, the following: paper, 

4 petroleum, chemical, pharmaceutical, food, and Coprocessing industries. 

5 The equations derived for the insertion probe are applicable to the cuvette mode of the 

6 ultrasound solute probe as shown in Figure 3. The transducers are identified as 22 and 24. In this 

7 case the cuvette 23 is housed in a controlled environment with a preset temperature. 

8 The clip-on mode of the ultrasound probe depicted in Figure 2 has a trough for the 

9 insertion of tubing 14. In this configuration, the transmission time from transducer 22 to 

10 transducer 24 includes the transmission time through the walls of the tube. However, the 

1 1 subtraction process employed in the derivation of Equation 1 1 also has the additional 

12 transmission time subtracted out so that Equation 1 1 remains valid for the probe of clip-on mode. 

1 3 In this embodiment, the clip 26 has a trough about 5 mm for the insertion of hemodialysis tubing 

14 14 whose outside diameter is about 6.2 mm. Ultrasound gel is used to facilitate the insertion of 

15 the tubing 14 into the clip and to provide an airtight contact between the tubing 14 and the 

16 transducers 22, 24, preventing problems and errors that can be caused by the reflection of 

1 7 ultrasound through air bubbles likely trapped between the transducers and tubing. 

1 8 Pipeline, tubing, or chromatography columns with a diameter other than 6.2 mm can also 

19 be used in the system and the dimensioning adjustments to the clip will be obvious to those 

20 skilled in the art. The resulting adjustment to convert phase shift to density or solute 

21 concentration will be incorporated through the parameters stored in a data file accompanying the 

22 probe. To gain a larger receiving signal, one can employ a lower ultrasound frequency such as 1 

23 Megahertz to power the emitting transducer. With the sampling frequency remaining at 65 
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1 MSPS, one will have 65 samples digitized over one period of oscillation while the total number 

2 of samples used in the determination of phase shift may remain in the range of 1000 to 2000. 

3 When the insertion or clip-on mode of the solute monitor is mounted onto the end of a 

4 column of chromatography process to assess solute concentration, we note that the transmission 

5 time includes that through the column wall, the porous beads in the column and the fluid. The 

6 presence of the beads will alter the relation between the fluid density and phase shift, which can 

7 be resolved by the on-line calibration procedure described later. 

8 In application, the ultrasound probe of insertion mode is inserted into a fluid processing 

9 system such as a pipeline or chemical processing system for continuous measurement of the 

10 phase shift of ultrasound transmission in the fluid. In addition, a test fluid with a density ptest at a 

1 1 volume AV is injected into the system upstream of the probe. Let the flow be Q and the volume 

12 of fluid situated between the injection and measurement site be Vi. According to the density 

13 dilution theory, one deduces from the principle of mass conservation the relation specified in 

14 Equation 17 

15 Q J(po - p)]dt = AV (po - ptest) Equation 17 

16 where p is the density being measured, po the steady state density before the injection, and the 

17 integration limit is over the time that the density is transiently deviated from the steady state 

1 8 density. Corrections can be made for the integration when the injected fluid re-circulates back 

19 through the probe. With most small injections, the density change from po is small. 

20 Equation 14 can be re-expressed as a linear relation between density change Ap and the 

21 phase change Acp with a calibration constant bi, i.e. 

22 Ap = bi A(p Equation 18 
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1 Its substitution into Eq. 17 yields Equation 19 for in situ^ on-line determination of the calibration 

2 constant 

3 bi = (AV/Q)([[?]]etest " [[?]]fio)/ J A(p dt Equation 1 9 

4 Using Equation 14, one can convert Equation 17 to Equation 20: 

5 Q = AV((ptest - <Po)/{/(9 -cpo)dt} Equation 20 

6 where [[?]]fio, [[?]]fitest, [[?]]fi, 9o, <ptest and (pare respectively the replacements of [[?]]fii, [[?]fi2, 

7 [[?]]£3, <P2 and 93 in Equation 14. Equation 20 can therefore be used for the calculation of the 

8 flow 

9 The consideration on mean transit time, volume and flow in the density dilution theory 

10 yields Equation 21 to relate the volume and flow: 

11 Q = Vj(9 - (po)dt /{J ((p - <po)tdt} Equation 21 

12 From Equation 20 and 21 , we can solve for the value of (ptest when the flow and volumes are 

1 3 known. The difference between the solute concentration in the flowing fluid and that in the test 

14 fluid can now be related to the difference between (ptest and (po- 

1 5 As another alternative to measure the flow, one can place two probes in two locations 

16 along the pipeline or chemical processing system. The flow can be calculated as: 

17 Q = V4-5J((p4 - (po)dt}/{ J(cp5 - (p4)tdt } Equation 22 

1 8 where V4.5 is the fluid volume of the pipeline or chemical system in between the two probes, 94 

19 is the phase shift measured by the upstream probe, and 95 that by the downstream probe. 

20 There are several working models of the ultrasound system to assess blood density. In 

21 one signal processing embodiment the transducer 22 of Fig. 1 and 5 is activated by any 

22 commercially available pulser/receiver, including a Panametrics 5072PR pulser/receiver, to emit 

23 ultrasound impulses of about 15 to 20 Megahertz at a rate of about 100 Hertz. Each pulse 
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1 contains about 4 to 6 oscillations. After its transmission through the flowing blood, the second 

2 transducer 1 8 receives the ultrasound pulse. The trigger signal from the pulser/receiver triggers a 

3 digital oscilloscope such as a LeCroy 9350AL oscilloscope or a Synatec A/D converter to sample 

4 the signal from the receiver 24 at 100 Megahertz to 1 Gigahertz. Some ten digitized samples of 

5 the pulse are compared by a computer for the determination of the phase shift. The sensitivity 

6 achievable with this process appears at best of the order nanosecond. 

7 In another embodiment, the signals from the receiving transducer and the function 

8 generator, which excites the emitting transducer, are processed by a phase lock amplifier such as 

9 Stanford Research System SR 844, which is powerful in processing signals with a frequency as 

10 high as 200 Megahertz. Many digitized samples from the continuous ultrasound are employed 

1 1 for the calculation. At the ultrasound frequency of interest here, the use of SR 844 provides 

1 2 sensitivity about one order of magnitude higher than the pulse system described previously. 

1 3 In our preferred embodiment, the ultrasound probe is used in conjunction with a novel 

14 ultrasound signal processing unit, a new computation scheme, and a two-fluid calibration 

15 procedure. This new scheme is an improved version of the system described in U.S. Patent 

16 Application No. 10/27 4 ,086 6,740,036 as we relax the requirement that the value of m in 

1 7 Equation 5 must be an integer multiple of 4. In the current invention, the digitized data of some 

18 100 oscillations (or 1000 to 2000 digitized samples) are employed. In comparison with the pulse 

19 procedure, the use of more samplings for signal processing enables us to gain higher sensitivity 

20 in phase shift determination. 

21 As in the SR 844 phase lock amplifier, the 90-degree shift of the excitation or receiving 

22 signal in Equation 4 is accomplished by electronic means. Our scheme achieves the shift by 

23 shifting the index of digitized receiving signals for multiplications and summations. Since the 
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1 shift in general may not be exactly 90-degree, the more general Equation 7 is deduced to carry 

2 out the computation of the phase shift between the emitting and receiving signal. Our test results 

3 indicate that our signal processing unit and the use of Equation 3, 4 and 6, even with a sampling 

4 rate of 65 MSPS (or about 15 ns a sample), can achieve a phase resolution of the order 0.03 ns, 

5 which is comparable to that via the SR844 amplifier. Only with this sensitivity, can the density 

6 and solute concentration expected to occur in industrial and clinical applications be measured. 

7 Because of the new algorithm and the simplicity in hardware design, our invention is in 

8 the position of being built as an IC chip. The end result of the use of IC chip will be a monitor 

9 with a much lower manufacturing cost than a system using SR 844 phase lock amplifier to 

10 process the signals. 

1 1 The solutes in the solution may exhibit different characteristics to absorbance or 

1 2 reflectance at various wavelength of the light. In a preferred embodiment, the ultrasound probe 

13 of the monitoring svstem works in conjunction with a light-emitting mechanism such an optical 

14 detector or a spectrophotometer to measure the absorbance and/or reflectance of fluid in the 

15 monitoring svstem. An optical detector preferablv comprises at least one light source and filter 

16 and at least two photomultipliers. The spectrophotometer preferablv comprises at least one light 

17 emitting diode and at least two photodiodes. C onventionallv, optical Conventionally, optical 

1 8 modularity modality requires the detector be calibrated with solutions of different solute 

1 9 concentration. However, as described later, one can employ the ultrasound probe and the on-line 

20 procedure to calibrate the optical modality. Let us illustrate this procedure with blood as the 

21 fluid and the optical modality being an IR detector. At an appropriate IR frequency, the detector 

22 has been used to assess hematocrit, the volumetric percentage of red blood cells in blood. 
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1 In clinical practice, the monitors are mounted onto the hemodialysis circuit, which 

2 withdraws blood from an artery and returns the blood after its passage of the hemodialysis 

3 machine back to a vein of the patient. A bolus of isotonic saline is injected into the circuit to 

4 flow through the hemodialysis machine and then the monitors. The sound velocity of saline is 

5 lower than that of blood. The passage of the saline after its mixing with the circulating blood 

6 will be recognized by the density monitor as a reduction in phase shift. Since there are no red 

7 blood cells in saline, we also expect to see a transient decrease in hematocrit, which will be 

8 picked up by the IR detector because of the resulting change in absorbance or reflectance 

9 characteristics of blood. Let us express the measurement of the optical detector as optical 

10 density (OD). When the blood is mixed with a saline of density ps and no red blood cells, the 

1 1 dilution of the density ([[?]]fib) and hematocrit (H) of blood follow Equation 21 . 

12 ([[?]]fib - ps) = [(pb,i - [[?]]fis)/Hi]H Equation 21 

13 where pb,i is the steady-state blood density and Hi is the steady-state hematocrit before the saline 

14 injection. 

15 By plotting the change in optical density AOD against that of density Apb as detected 

16 through the ultrasound probe, we can obtain the slope b2 in the linear relation of Equation 22: 

17 Apb= b2 AOD Equation 22 

18 Its combination with Equation 21 for saline dilution yields Equation 23: 

19 AH/Hi = [b2/(pb,i - ps )]AOD Equation 23 

20 In practice, the optical detector is located, for example, downstream of the ultrasound probe. To 

21 carry out the data analysis through Equation 22, we will adjust the optical signal by a time for 

22 which the linear fit between the optical density and density has the correlation coefficient closest 

23 to unity. Then the slope of this linear fit is taken as b2. Once the slope is measured with this on- 
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1 line procedure. Equation 23 is the calibration equation to convert the change in optical density to 

2 the hematocrit ratio (AH/Hi). 

3 In industrial applications, the injectate may contain a number of solutes of interest to the 

4 chemical and biological process. Each solute may have different absorbance or reflectance 

5 characteristics. Thus the spectrum of the optical detector will be set up to differentiate the 

6 absorbance or reflectance of the solutes. Because of difference in molecular sizes or activities, 

7 the passage of these solutes through the chemical system, biological reactor or chromatography 

8 column may occur at various times. Their presence in the flowing fluid will alter the phase shift 

9 and light absorbance as the fluid passing through the monitor. A procedure similar to hematocrit 

10 calibration can be applied as an on-line means to calibrate the optical detector in industrial 

1 1 setting. Conversely, the specificity of the optical detector in solute identification allows us to 

12 pinpoint which passage detected by the ultrasound probe is associated with which specific solute. 

13 The passage of solute as detected by the ultrasound probe can now be used to activate a 

14 fractional collector to collect the solution containing most of the specific solute. This on-line 

15 control will reduce the collection of other solutes coming through the probe at other times and to 

16 improve the purity of the specific solute in the collection. This characteristic to identify the 

1 7 presence of solutes through the use of density/solute monitor can also be used to improve the 

1 8 collection of solution containing no solutes for reuse by the chromatography process. 

1 9 Y **^y is absorbed by the fluid over its passage. The attenuation of y ray, a physical 

20 characteristic of the fluid, can be associated with and thus be used to determine the density of the 

21 fluid. The combination of this attenuation with the ultrasound characteristics may further 

22 enhance our ability to differentiate the kind of fluid flowing pass the density/solute monitor. 
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1 Infusion of hypertonic saline has been used cHnically to extract fluid from the tissue in 

2 human body to the circulation. The extracted fluid has a density lower than the blood. Thus its 

3 mixing with blood will lower the density of blood. Consequently, multi-modality monitor on 

4 conductivity and phase shift may allow us to assess the process of fluid extraction from the 

5 tissue. 

6 Tracers have been used to monitor dynamic events occurring in a chemical system, 

7 pipeline, or the human circulation system. Some tracers can be detected by magnetic resonance 

8 imaging (MRI) system or CT scan. If the tracers are in the form of vesicles containing a fluid or 

9 other medium that its density is different from that of the flowing fluid, then the density or 

10 compressibility of fluid may be altered by the presence of tracers and be detected by the density 

1 1 monitor. The tracer can be a substance tagged with radioactive element or dye. Radioactivity 

1 2 counter can detect the former and optical detector set at the frequency most sensitive to the dye 

13 can detect the latter. The use of multi-modality detection systems and tracers may further 

14 improve the sensitivity and specificity of the density/solute monitor to better track the movement 

15 of solutes through chemical processing system or human circulation. 
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1 Abstract 

2 A density/solute monitor having at least one ultrasound probe, a signal processing unit, 

3 and a computing mechanism, and process for using the same, to measure phase shift between 

4 emitting and receiving ultrasound, sound velocity, compressibility, density, and solute 

5 concentration of fluid flowing through a fluid processing system. The ultrasound probe emits 

6 and receives ultrasound waves through the fluid and the signal-processing unit and computing 

7 mechanism process the ultrasound waves to determine phase and time shift. The computing 

8 mechanism converts phase shift to density, compressibility, and solute concentration 

9 measurements of the fluid. Calibrating fluids calibrate the detected phase shift in terms of sound 

10 velocity in the factory. Measurements provide information about passage of solutes and flow to 

1 1 achieve better solute collection efficiency, solution purity, and control of fluid processing 

12 systems. The density/solute monitor can include other detection modalities such as an optical 

1 3 probe, making concentration measurements of the density/solute monitor more specific to 

14 selected solutes. 
15 
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